Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



Chandra observation of an unusually long and intense X-ray flare 

from a young solar-like star in M 78 

N. Grosso 1 ' 2 , T. Montmerle 2 , E.D. Feigelson 3 , and T.G. Forbes 4 



1 Max-Planck-Institut fur extraterrestrische Physik, P.O. Box 1312, D-85741 Garching bei Munchen, Germany 

2 Laboratoire d'Astrophysique de Grenoble, Universite Joseph-Fourier, F-38041 Grenoble cedex 9, France 

3 Department of Astronomy and Astrophysics, Pennsylvania State University, 525 Davey Laboratory, University 
Park, PA 16802, USA 

4 Institute for the Study of Earth, Oceans, and Space, 39 College Road, University of New Hampshire, Durham, 
NH 03824, USA 

Received 7 July 2003 / Accepted 29 January 2004 

Abstract. LkHa312 has been observed serendipitously with the ACIS-I detector on board the Chandra X-ray 
Observatory with 26 h continuous exposure. This Ha emission line star belongs to M 78 (NGC 2068) , one of the 
star-forming regions of the Orion B giant molecular cloud at a distance of 400 pc. From the optical and the near- 
infrared (NIR) data, we show that LkHa312 is a pre-main sequence (PMS) low-mass star with a weak NIR excess. 
This genuine T Tauri star displayed an X-ray flare with an unusual long rise phase (~8h). The X-ray emission 
was nearly constant during the first 18 h of the observation, and then increased by a factor of 13 during a fast 
rise phase (~2h), and reached a factor of 16 above the quiescent X-ray level at the end of a gradual phase (~6h) 
showing a slower rise. To our knowledge this flare, with ~0.4-~0.5 cts s -1 , has the highest count rate observed 
so far with Chandra from a PMS low-mass star. By chance, the source position, 8.2' off-axis, protected this 
observation from pile-up. We make a spectral analysis of the X-ray emission versus time, showing that the plasma 
temperature of the quiescent phase and the flare peak reaches 29 MK and 88 MK, respectively. The quiescent and 
flare luminosities in the energy range 0.5-8 keV corrected from absorption (Ah ~ 1.7 10 21 cm -2 ) are 6 10 30 ergs -1 
and ~10 32 ergs -1 , respectively. The ratio of the quiescent X-ray luminosity on the LkHa 312 bolometric luminosity 
is very high with log(Lx/£boi) = —2.9, implying that the corona of LkHa 312 reached the 'saturation' level. The 
X-ray luminosity of the flare peak reaches ~2% of the stellar bolometric luminosity. The different phases of this 
flare are finally discussed in the framework of solar flares, which leads to the magnetic loop height from 3.1 10 10 
to lO^cm (0.2-0.5 R*, i.e., 0.5-1.3 R©). 

Key words. Open clusters and associations : individual : M 78, NGC 2068 - X-rays : stars - Stars : individual : 
LkHa312 - Stars : flare - Stars : pre-main sequence - Infrared : stars 



1. Introduction 

Since the first X-ray observations of star-forming re- 
gions (SFR) with the Einstein Observatory, young 
PMS low-mass stars, T Tauri stars (TTS), are known 
to be variable in X-rays ( |Feigelson fc Decampli 1 981 
IMontmerle et al. 1983J) . TTS, as active stars, dis- 
play X-ray flares triggered by magnetic reconnection 
events occurring in their stellar coronae (see reviews by 
|Feigelson fc Montmerle T999| and lFavata fc Micela 2003|l . 
The Chandra X-ray Observatory, thanks to its highly 
elliptical orbit that permits continuous observation over 
many hours, has collected a real zoo of X-ray flares from 
TTS in several SFR (e.g., IPreibisch fc' Zinnecker 20021 



Feigelso n et al. 2002a| llmanishi et al. 2003). The X-ray 
flares of TTS are generally impulsive, displaying a fast rise 
phase (~2h) corresponding to a heating phase, followed 
by an exponential decay corresponding to a cooling phase 
(e.g., fmanishi et al. 2003). Only a few long duration 
flares have been observed from TTS with the previous 
generation of X -ray satellites, ROSA T and ASCA (SR13, 
ICasanova l"99l object V 7 73Tau, ISkinner et al. 19971 
P1724, IStelzer et al. 1999|l . |Stelzer et al. (1999)| showed 
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that these light curves can be reproduced by the rota- 
tional modulation of the exponential decay of a cooling 
flare. Feigelson et al. (2002a)| investigating the variabil- 
ity of young solar-like stars (M*=0.7-1.4M Q ) in the 
Orion Nebula Cloud, find several likely long-duration 
flares but the low count rates of these events prevent a 
time-dependent spectroscopy analysis. 
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Fig. 1. Chandra observation of LkHa 312 on October 18, 2000. The background image is the ACIS-I observation in 
the energy range 0.5-8 keV of LkHa 312, located at ~ 8' off-axis (the pixel size is 0'.'25). The color scale is stretched 
logarithmically to enhance the point spread function (PSF) features, note the 2" asymmetric cusp surrounded by a 
10" elliptical halo. The contour map displays arbitrary intensities from the 2MASS J-band image using linear scale. 
The positions of X-ray sources detected above the threshold of 3 10 -7 are marked with 'x'. Their infrared positions 
in the 2MASS survey are marked by '+'. 



M78 (NGC2068) is a reflection nebula il- 
luminated by a B1.5 V star (HD 38563 North; 
IMannion fc Scarrott 19 84). and located in the northern 
part of the closest giant molecular cloud, Orion B (L 1630), 
at a distance of ~400pc | |Anthony-Twarog 1982| ). The 
star-forming region M 78 has been observed with the 
ACIS-I detector aboard Chandra with 26 h exposure on 
October 18, 2000 (sequence number 200100). We report 
here the study of the brightest X-ray source detected 
during this observation, LkHa 312, located ~ 10' at 
the South-East of the optical emission nebula M78. 
This source displays an intense X-ray flare with an 
unusually long rise phase. We describe the ACIS-I data 
reduction and the X-ray detection of LkHa 312 in |J21 
The evolutionary status of LkHa 312 is determined from 
the optical and the NIR data in <J|1 The Chandra light 
curve of the flare is presented in 21 The time-dependent 
spectra is investigated in SjJJ Finally, the different phases 
of this flare are discussed in the framework of solar flares 
in Concluding remarks are presented in 

2. Chandra ACIS-I observation 

The telemetry format of this observation is 'Timed 
Exposure' with 3.2 s for the frame time, combined with the 
ACIS mode 'Faint'. The aimpoint of the ACIS-I detector is 
aj2ooo = 05 h 46 m 44?2, S J20 oo = +00°03'46'/0 (I = 205.3°, 
b = —14.3°), centered on the optical emission nebula M 78. 

Following Feigelson et al. (2002b) we start our data 
analysis from the level 1 processed event list pro- 
vided by the Chandra X-ray Center (CXC) with the 



Table 1. Chandra X-ray sources with 2MASS counterpart 
in the vicinity of LkHa 312 and their X-ray count rates on 
October 18, 2000. 



CXOU 0547 


Count rate a 


2MASS J0547 


r b 


LkHa 




[ctsks -1 ] 




["] 




06.9+000047 


6.8+0.3 


0696+0000477 


0.3 


309 


08.7+000013 


5.7+0.3 


0868+0000140 


0.5 




14.0+000015 


149.8+1.3 


1408+0000157 


0.1 


312 



a In the energy range 0.5-8 keV, from wavdetect taking 

93859 s for the exposure. 
b Offsets between X-ray and NIR positions. 



pipeline processing version R4CU5UPD11. We first use 
the tool acis_process_events of the Chandra Interactive 
Analysis of Observations (CIA0) package (version 2.2.1 
and 2.3 1 ) to remove the unnecessary ±0.25" randomiza- 
tion introduced by the pipeline processing on each event 
position, and we select events belonging to Good Time 
Intervals (GTI) provided by the CXC. The observation 
time corrected from deadtime (after each CCD expo- 
sure 0.04 s are lost to read the CCD chip) is 93859 s, 
i.e., more than 26 h. Second we correct the energy and 
grade of each event using the Penn State's Charge 
Transfer Inefficiency (CTI) corrector package (version 
1.41 2 ; |Townsley et al. 2000] |Townsley et al. 2002| ), which 



l^p^/as^^arvard^du/c^o| 
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Table 2. ROSAT/HRI observation of LkHa312 from 
September 5 to October 5, 1997. 



ROR 


Exp. 


RX J0547 


£ a 


CR b 


2025 


M 






[cts ks _1 ] 


21h-l 


26.7 


14.1+000011 


213 


8.3±0.7 



C is the likelihood of existence. 

The count rate is for the ROSAT energy band, 0.1-2.4 keV. 



models the CTI characteristics of each CCD chip. The 
most fundamental result of this CTI correction is the im- 
provement in spectral resolution. Third we clean the data 
by selecting only events with ASCA grade (0, 2, 3, 4, and 
6) and energy in the range 0.5-8 keV after CTI correc- 
tion. Bad pixels and hot columns are rejected by selecting 
only events with status bits 1-15 and 20-32 equal to zero. 
Remaining events with non zero status bits are events 
which were identified by the pipeline processing as pro- 
duced by cosmic ray afterglows. However few of them are 
obvious events from bright X-ray sources misclassified by 
the pipeline. Hence we reject these events only during X- 
ray source detection to avoid spurious X-ray sources, but 
we use them for source spectral analysis. Finally to im- 
prove the absolute event positional accuracy, we correct 
the aspect offset using the fix offsets thread 3 . 

Fig- ID presents an enlargement of the resulting X- 
ray image at an off- axis angle 9 ~ 8'. Three X- 
ray sources counterparts of optical emission-line stars 
QHerbig fc Kuhi 1963) and/or NIR 2MASS sources (all- 
skv data release: ICutri et al. 2003}) are visible. To ob- 
tain accurate X-ray positions and to detect weaker X-ray 
sources in this area, we use the tool wavdetect of the CIA0 
package, with the threshold probability of 3 10~ 7 , and 11 
wavelet scales ranging from 0.5" to 16" in power of a/2. 
Table ^ gives the list of the Chandra X-ray sources in 
Fig. [T] identified with 2MASS counterparts. We find that 
the agreement between the X-ray and the NIR positions 
is better than / .'5. No X-ray source is detected at the 
position of LkHa313, but 4" away we find the remark- 
able X-ray detection of LkHa312 with more than 14,000 
counts, i.e., an average count rate of ~0.15ctss _1 . 

We find a previous X-ray detection of this source made 
from September 5 to October 5, 1997, with the High 
Resolution Imager (HRI) on board ROSAT in the energy 
range 0. 1-2.4 keV. Table [2] gives the results of our analy- 
sis of the archived ROSAT Observation Request (ROR) 
with the EXSAS 4 package l|Zimmermann et al. 1997jl . The 
total exposure of 26.7 ks is the sum of several ~ 20min 
segment observations spread over one month. We extract 
source events from the 95% encircled energy radius and 
produce a background subtracted light curve (corrected 
from the PSF fraction) with one bintime per segment 
(Fig. [5J. As the HRI has no spectral resolution, we es- 
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Fig. 2. ROSAT/HRI X-ray light curve of LkHa312 ob- 
tained from September 5 to October 5, 1997. Crosses show 
the X-ray light curve and one-sigma error bars in the 
energy range 0. 1-2.4 keV for each observational segment 
(^20min). Background count rates are shown as upper 
limits when the background subtracted count rate is neg- 
ative. 

timate the ROSAT/HRI count rate of the quiescent level 
with W3PIMMS 5 using the Chandra light curve and spec- 
troscopy results (see low inland |J5J, which give a Chandra 
quiescent count rate of 0.035 cts s -1 , and a Chandra spec- 
trum model with N H = 1.73 10 21 cm" 2 , Z ~ 0.2 Z Q , 
< kT >=2.5keV. We find 0.004 cts s _1 , which is consis- 
tent with the low level detected visible in the light curve. 
This source is clearly variable, however the poor sampling 
of the light curve precludes obtaining its detailed shape. 

3. Evolutionary status of LkHa 312 

LkHa312 was originally selected due to its (weak) 
Ha line emission in slitless grating survey of M78 
( |Herbig fc Kuhi 1963| ) . An optical spectrum with 7 A reso- 
lution has been obtained in 1977 by Cohen fc Kuhi (1979)1 
leading to its M0 spectral type classification, and its Ha 
equivalent width, EW(Ha)=8.2A (see Table EJ|. We note 
that this star is not classified as a TTS in the catalog 
of |Herbig fc Bell (1988)1 probably due to its low EW(Ha) 
value. |Martm (1997)| considers that classical TTS (CTTS) 
should have EW(Hq) larger than expected from the av- 
erage chromospheric activity for their spectral type, and 
adopts EW(Ha)=10A for early-M stars as a lower limit, 
instead of the classical EW(Ha)=5A limit. On the ba- 
sis only of its EW(Ha), LkHa312 could then be either 
an non-accreting PMS star of M 78 or a foreground main- 
sequence dwarf with Ha emission, a dMe star. Of course 
the detection of the Li I 6707 absorption line would defini- 
tively demonstrate the youthfulness of this star, but this 
observation is not available. 



http : //exc . harvard . edu/cal/ ASPECT/f ix_of f set/f ix_of f 
|h^p^77wave^crayT mpe .mpg.de/ex sasj 



Recently Preibisch & Zinnecker (2002) have shown in 
the young star cluster IC 348 that Ha emission is a poor 
indicator of circumstellar material around TTS. More pre- 
cisely, an IR excess in the L-band or beyond may be 

nrpspnt at the same time as weak Ha- emission, indicat- 
set . cgi 

5 http : //asc .harvard. edu/toolkit/pimms . jsp 



4 



Grosso et al.: An unusually long and intense X-ray flare from a young solar- like star in M78 



Table 3. Optical and available near-infrared data for the emission line stars discussed in this article. See corresponding 
J — H, H — Ks diagram in Fig. [SJ spectral energy distributions in Fig. and H.-R. diagram in Fig. 
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This work 
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T e ff 
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J a 




K s * 
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T b 


# 
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[A] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[Lq] 


309* 


K7 


4000 


51.2 


15.1 


11.35±0.04 


10.57±0.04 


10.22±0.04 


0.78±0.08 


0.35±0.08 


2.0±0.5 


1.6±0.2 


312 


MO 


3917 


8.2 


15.1 


11.32±0.05 


10.52±0.06 


10.12±0.06 


0.80±0.11 


0.40±0.12 


1.8±0.5 


1.5±0.3 


313* 


MO. 5 


3802 


19.6 


15.9 


10.75±0.09 


9.79±0.09 


9.23±0.07 


0.95±0.18 


0.57±0.16 


3.5±0.5 


3.9±0.5 



a The errors are the total magnitude uncertainties given at the 90% confidence level. 

b From fits of Fig. 0] and taking d = 400 pc. 

* Classified as classical T Tauri star in Herbig & Bell (1988). 




Fig. 3. J — H, H — K s diagram of LkHa 312 (open dot). 
Magnitude error bars are the total magnitude uncertain- 
ties at the 90% confidence level. The intrinsic colors of 
giants (dotted line) and A0-M6 dwarfs (continuous line) 
from Bessel & Brett (1988), adapted for the 2MASS pho- 
tometric system using the 2MASS color transformations 
(Carpenter 2001; Cutri et al. 2003), are plotted for com- 
parison. The extinction law of Cohen et al. (1981), also 
adapted for the 2MASS photometric system, is drawn for 
giants and M9 stars. The arrow indicates a visual extinc- 
tion of 1 mag. Black dots mark the position of known clas- 
sical T Tauri stars in the vicinity of LkHa 312 (Herbig & 
Bell 1988). 



ing that such objects are in fact surrounded by a hollow 
circumstellar disk but currently are accretion quiet. A sim- 



ilar conclusion was drawn by Andre & Montmerlc (1994) 
based on a sample of p Ophiuchi dark cloud's TTS. 
Preibisch & Zinnecker (2002) use K — L to trace the in- 



frared excess, but L-band photometry is not available for 
LkHa 312. We have found in the ISO archives 6 , 6.7/im 
and 14.3 /im ISOCAM detections of the TTS LkHa 313, 
but the poor angular resolution of ISOCAM (~6") can- 
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Fig. 4. Spectral energy distributions. The optical spec- 
tra and the NIR photometry come from Cohen & Kuhi 
(1979) and 2MASS, respectively. Dashed line shows the 
absorbed black body spectrum having the effective tem- 
perature found by (Cohen & Kuhi 1979), coinciding with 
the J-band photometry, and consistent with the optical 
spectra. Upper and lower dotted lines show the same 
black body spectrum for minus and plus 0.5 mag extinc- 
tion, respectively. We give the resulting bolometric lumi- 
nosities, including typical ranges for a typical extinction 
error of 0.5 mag. 

not resolve the emission of LkHa 312 only 4" away. Fig- EI 
displays in a J — H, H — Ks diagram the position of two 
known TTS of M 78 selected due to their strong Ha line 
emission in the vicinity of LkHa 312 (see Table EJ. Both 
TTS have a visual extinction, Ay, greater than 1 mag. The 
NIR colors of LkHa 312 are very similar to those of the 
accreting TTS LkHa 309, and can be retrieved from the 
intrinsic color of a M0 star by applying ~ 1.8 mag of vi- 
sual extinction combined with a small NIR excess, which 
is maybe produced by a circumstellar disk with an inner 
hole IjAndre fc Montmerle 1994|l . This extinction is con- 
sistent with LkHa 312 being a young star of M78 rather 
than a foreground main-sequence dwarf. 
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We note that the extinction derived from the color- 
color diagram for LkHa 312 is a factor of ~4 larger 
than the one found by Cohen fc Kuhi (1979)| from op- 
tical spectroscopy. The extinction values obtained by 
Cohen fc Kuhi (1979)| for LkHa 309 and LkHa 313 are 
also smaller than the one measured in our color-color dia- 
gram. Moreover the bolometric luminosities. Lboli given by 
Cohen fc Kuhi (1979)| for these objects "represents solely 
the optical data (extrapolated [with a black body] to 
infinite wavelength from 0.67 fim)" (see their Table 5, 
note of columm 17). To investigate this extinction dis- 
crepancy and to udpate the bolometric luminosities, we 
build spectral energy distributions (SED; Fig. 0J. The 
optical spectra were copied directly from the Fig. 23 



10.0 



of |Cohen fc Kuhi (1979)] The 2MASS photometry was 
first transformed to the UKIRT system (Carpenter 2001; 
Cutri et al. 2003), and then converted to flux densities 
( Co x 2000jl . We compute the black body spectrum having 



the effective temperature found by Cohen & Kuhi (1979) 



we apply using |Fitzpatrick (1999) the extinction given by 
Cohen fc Kuhi (1979)1 and" 

we scale the absorbed black 
body spectrum for coincidence with the 0.67 /im photom- 
etry. The obtained SED are far below the NIR data for 
all the stars of Table hence both the extinction and 
the bolometric luminosity of Cohen & Kuhi (1979) are un- 



derestim ated, and need to be u dpated. Taking the value 
found by Cohen fc Kuhi (1979) for the effective tempera- 
ture, we adjust the extinction to obtain an absorbed black 
body SED, coinciding with the 2MASS photometry in the 
J-band (the contribution of infrared-excess, if any, is neg- 
ligible at this wavelength) , and consistent with the optical 
spectrum. Fig. 01 displays our best absorbed black body, 
with the resulting Lboi estimate. We consider a typical 
Ay error of 0.5 mag. We list in Table |3 our values of Ay 
and Lboi for LkHa 309, LkHa 312, and LkHa 313. The es- 
timated extinctions are now consistent with the ones de- 
rived from the color-color diagram. 

Fig.JSJshows the H.-R. diagram with the PMS tracks of 
Baraffe et al. (1998) fitting the Sun (i.e., having the mix- 
ing length parameter, / m i x , equal to 1.9 times the pres- 
sure scale height, H p ), with the position of LkHa 312 
assuming that LkHa 312 and the other stars of Table 
are at the distance of M 78 (d ~ 400 pc) . Considering 
LkHa 312 as a main sequence star would put it on the 
120 Myr isochrone revising down its distance to only 
d ~ 90 pc. This would locate LkHa312 in the Local 
Bubble where N K ~ 510 19 cm" 2 IjSfeir et al. 1999jl . This 
small distance would not be consistent with the extinc- 
tion of this object, Ay ~ 1.8 mag, which implies Nn ~ 
3.210 21 cm -2 ( using, e.g., the Nn/Ay conversion fac- 
tor of Pred ehl fc Schmitt 1995|l . Therefore we consider 
LkHa 312 as a PMS star of M 78. 



From the |Baraffe et al. (19 98) PMS tracks we find that 
LkHa 312 is located above the lMyr isochrone, among 
known CTTS of M 78, with a stellar mass of 0.7-0.75 M Q . 
We also note that if LkHa 313 were a close binary sys- 
tem the luminosity of the individual components would 
be decreased by a factor of ~2 (see Fig. 0). In this case 
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Fig. 5. H.-R. diagram of LkHa 312 (open dot). Dotted 
lines show Baraffe et al. (1998) PMS tracks (with the mix- 
ing length required to fit the Sun, / m i X = 1.9 H p , with 
Hp the pressure scale height). Dashed lines show the cor- 
responding isochrones. Black dots mark the position of 
known classical T Tauri stars in the vicinity of LkHa 312 
(Herbig & Bell 1988). The dashed circle marks the posi- 
tion of LkHa 313 in case of binarity. Error bar shows the 
typical error on the effective temperature (Cohen & Kuhi 
1979), and on the bolometric luminosity (see Fig.^J. 



LkHa 309, LkHa 312, and LkHa 313 would all be coeval, 
and LkHa 312 and LkHa 313 would form a gravitationally 
bound system since their separation is only 4" on the sky. 

Our final conclusion on the evolutionary status of 
LkHa 312 is that it is a genuine TTS, with maybe a transi- 
tion circumstellar disk, belonging to the M 78 star-forming 
region. Combining T c g = 3917 K and Lboi = 1-5 Lq leads 
to the stellar radius i?* = 2.6 R Q . 



4. Variability study 

4.1. Chandra X-ray light curve 

We extract source events from the 99% encircled energy 
radius, R(99%EE). We use the empirical fit given by 
Feigelson et al. (2002b)| R(99%EE) = 8 + 0.20 (with 
9 the off-axis position), which leads for 8 = 8.2' to 
R(99%EE) = 9.6". The background is taken from a larger 
area free of sources. We use the CIA0 tool lightcurve, 
which takes into account Good Time Intervals, to pro- 
duce background subtracted light curves (corrected from 
the PSF fraction) with 20min bintime, and energy filter- 
ing. 

Fig.Elshows the light curve of LkHa 312 in the energy 
range 0.5-8 keV. The X-ray emission was nearly constant 
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Fig. 6. Chandra background subtracted X-ray light curves of LkHa 312. The time binning is 20 min. Crosses show 
the X-ray light curve and one-sigma error bars in the energy range 0.5-8 keV. The stripes show the X-ray light curves 
including onc-sigma errors corresponding to the energy ranges 0.5-1.6keV (soft band), 1.6-2.6keV (hi band), and 
2.6-8 keV (h 2 band). Vertical dashed lines and horizontal arrows define the time intervals, Ati, used for the time- 
dependent spectroscopy. The quiescent, fast rise, and gradual phases correspond to the time intervals Ato, Ati-At2, 
and At3-Atg, respectively. 
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Fig. 7. Left panel : hardness ratios of LkHa 312 versus time. The top and bottom panels show respectively the hardness 
ratio HR1 and HR2, computed from the counts collected in three energy bands (soft, hi, and h 2 , see Fig. EJ), with 
1.64cr error bars. The grey stripe indicates the average hardness ratio of the quiescent phase. Right panel : hardness 
ratios of LkHa 312 compared to spectral models. Iso-energy and iso-absorption lines correspond to absorbed optically 
thin one-temperature plasma spectra model with 0.2 solar abundances, and show the correspondence between the 
hardness ratio pair (HR1, HR2) and the absorption/temperature pair (Nh, T). The grey area shows the average of 
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during the first 18 h of the observation, defined as the time 
interval Ato, corresponding to the quiescent phase. The 
count rate was then multiplied by a factor of 13 during 
a fast rise phase (~2h) corresponding to Ati-At2, and 
reached a factor of 16 above the quiescent X-ray level at 
the end of a gradual phase (~6h) showing a slower rise 
and corresponding to Ata-Atg. The duration of the fast 
rise phase is similar to the one observed in other TTS (e.g., 
Ilmanishi et al. 20f)3J) . however the gradual phase is really 
unusual. The light curve shape is reminiscent of one of 
the egress phase of an eclipse; however, as we will see be- 
low, this event was triggered by plasma heating. Therefore 
LkHa312 displayed a long duration X-ray flare, in con- 
trast to an impulsive flare, during the phase Ati-Atg. 
Some roughly similar events were already observed by 
Chandra in the Orion Nebula Cluster (see for example, 
JW567 and JW738 in Fig. 3 of |Feigelson et al. 2002aD , 
or in IC348 (CXOPZ J034 416.4+320955 in Fig. 3 of 
IPreibisch <k Z innccke r~2"002Jl but the X-ray counts were 
lower, and the datastream was shorter. 

The flare count rate was very high, ~0.4-~0.5ctss -1 , 
allowing time-dependent spectroscopy (see jjjjj. To our 
knowledge this event has the highest count rate observed 
so far from a PMS low-mass star with Chandra. Such a 
high count rate implies a pile-up fraction of ~20% for an 
on-axis source, however for 9 = 8.2' the PSF spreading 
reduces the pile-up fraction to a value lower than 1% (see 
The Chandra Proposers' Observatory Guide 7 ), and hence 
we do not need to worry about this effect in our spectral 
studies. 

4.2. Hardness ratios 

Before applying automatic fitting procedures to spectra 
in §5\ we investigate count ratios ('hardness ratio' equiv- 
alent to 'X-ray colors') variation, to assess absorption 
and/or temperature changes during the flare. We split 
the 0.5-8.0 keV energy band into three parts : the soft 
band, the hi, and h 2 hard bands, corresponding to the 
energy ranges 0.5-1.6keV, 1.6-2.6keV, and 2.6-8 keV, re- 
spectively. Fig. [5] shows the light curves of these three dif- 
ferent bands. The global shape of these light curves is very 
similar to the full energy band light curve. We compute the 
following hardness ratios : HR1 = (hi +h% — soft) /(hi + 
h 2 + soft), and HR2 = (h 2 — hi) /(hi + h 2 ), with asso- 
ciated errors computed using Gaussian error propagation. 
The left panel of Fig. shows the flare hardness ratios 
versus time. They display higher values during the flare 
phase than during the quiescent phase. 

We use rmf and arf files of the quiescent phase of 
LkHa312 (see method to compute these files below 
to simulate with XSPEC synthetic Chandra spectra of ab- 
sorbed (wabs model with IMorrison fc M cCammon 1983 
cross-sections) optically thin one-temperature plasma 
(mekal model; IKaastra et al. 1996jl with Z = 0.2 Z Q 
(found from the spectral fitting of the quiescent phase, 



sec below *15. 21 Tabic l4l. for a grid of absorption (N-r) and 
temperature (T) values, and we compute their correspond- 
ing hardness ratios. The right panel of Fig.|7|shows the re- 
sulting spectral model grid in the hardness-space with iso- 
energy and iso-absorption lines. Indeed the energy bands 
were designed to have the best separation of these iso- 
energy and iso-absorption lines, which allow a direct cor- 
respondence between the hardness ratio pair (HR1, HR2) 
and the absorption/temperature pair (JVh, T). 

The quiescent phase corresponds in this one- 
temperature plasma model to log-/Vn ~ 20.7, and kT ~ 
2.7 keV. However as will be shown below in |S]the quies- 
cent phase is in fact better modeled by a two-temperature 
plasma. The plasma temperature found with the hardness 
ratio method is consistent with the average of the plasma 
temperatures weighted over the emission measures, but 
the absorption estimate is here clearly underestimated, 
because as it will see below in (Table the spectral 
fitting leads to log-ZVu = 21.2 with a 90% confidence in- 
terval logJV H =21. 1-21.4. We add in Fig. |7|our best esti- 
mate of the (iVn, T) pair for the quiescent phase from our 
spectral fitting. We can only conclude with these hardness 
ratios that during the flare the temperature of the plasma 
increases to ^7keV on average. There are no variations of 
the column density. 

5. Time-dependent spectroscopy 

5.1. Methodology 

The spectral analysis is performed with XSPEC (version 
11.2 s ). Plasma models are convolved with the response 
matrix file (rmf) describing the spectral resolution of the 
detector as a function of energy, and an auxiliary response 
file (arf) describing the telescope and ACIS detector ef- 
fective area as a function of energy and location in the 
detector. Reduction in exposure times due to telescope 
vignetting, and possible bad CCD columns moving on 
the source with satellite aspect dithering, is included in 
the arf file. We use the rmf files and the quantum ef- 
ficiency uniformity (QEU) files provided with the CTI 
corrector package. These QEU files are used by mkarf to 
generate arf file consistent with our CTI corrected events. 
We note that for LkHa312 (II chip, row 749) far away 
from the readout node (II chip, row 1), the CTI correc- 
tion improves the energy resolution from FWHM~130eV 
to 110 eV in the soft part of the spectrum, and from 
FHWM~2 60eV to 190 eV in th e hard part of the spec- 
trum (see |Townsley et al. 2002| ). We apply the acisabs 
model to the arf file to take into account the contami- 
nation of the ACIS optical blocking filter. For our obser- 
vation 453 days after launch, this correction is negligible 
above ~1 keV, but leads to a decrease of 20% of the effec- 
tive area at 0.5 keV, which is not negligible to determine an 
accurate value of hydrogen column density on the line of 
sight. Best- fit model parameters are found by comparing 



7 http : //asc .harvard. edu/udocs/docs/POG/MPOG/index. html 8 http : //lheawww.gsfc .nasa.gov/users/kaa/home .html 
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Table 4. Best fit parameters of the time-dependent spectroscopy with their errors at the 90% confidence levels 
(Ax 2 =2.71; corresponding to a = 1.64 for Gaussian statistics). r\ = log (£x,intr/£boi)- 
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models with the pulse height distribution of the extracted 
events using x 2 minimization. The source spectrum is re- 
binned with grppha in the f tools package to a minimum 
of 20 source counts per spectral bin. The model is im- 
proved until obtaining xt ~ 1 an d 2 — 5%, with xt the 
reduced x 2 (i- e -; X 2 over the degree of freedom, v), and 
Q the probability that one would observe xt, or a, larger 
value, if the assumed model is true, and the best-fit model 
parameters are the true parameter values. The parameter 
error at the 90% confidence level (A% 2 =2.71; correspond- 
ing to a — 1.64 for Gaussian statistics) is computed with 
the commands error and steppar. 

5.2. Multi-temperature plasma model 

We first model the plasma in XSPEC with an optically 
thin single-temperature plasma model (mekal model; 
IKaastra et al. 1996)) . combined with photoelectric absorp- 
tion using Morrison & McCammon (1983) cross-sections 



(wabs model). We are led to increase the number of plasma 
component to improve the goodness-of-fit. 

5.2.1. Quiescent emission and flare spectra 

Fig. |H1 shows the quiescent and the flare spectra for com- 
parison. Both phases require more than a one-temperature 
plasma model to be fitted, and we use a two-temperature 
model. The resulting fit parameters are listed in Table 
The quiescent spectrum is featureless, no emission lines 
are visible. Approximatively 70% of the X-ray emission 
comes from a high temperature plasma at energy kT = 
3.1 keV. The plasma abundance is Z ~ 0.2 Zq (this 
is the value that was used previously for hardness ra- 
tio study in H4.2|) . The average plasma energy, weighted 
over the emission measure, is < kT >=2.5keV, cor- 
responding to an average temperature of 29 MK, sim- 



ilar to the temperature found in the previous section 
with the hardness ratio method. The absorption is low 
with JVh = 1.3-2.3 10 21 cm -2 corresponding to Ay -0.7- 
1.3 mag IjPredehl fc Schmitt 1995) . 

The intrinsic luminosity in the energy range 0.5-8 keV 
is ix,intr = GlO^ergs -1 . From the stellar radius i?* = 
2.6 Rq inferred from the H.-R. diagram (see above Q, we 
find a high X-ray stellar flux, Fx = 1-4 10 7 erg cm -2 s _1 . 
To quantify the X-ray activity we use the logarithm of the 
ratio of the X-ray luminosity corrected from absorption 
on the bolometric luminosity, 77 = log(Fx,mtr/Fboi)- (This 
indicator is independent of the source distance.) The qui- 
escent phase displays 77 = —2.9. This value is very high 
compared with the average value for TTS, r\ ~ —4, and 
implies that the corona of LkHa312 reached the satura- 
tion level. 

The large number of events collected during the flare 
phase (about 12,000) unveils in the flare spectral resid- 
ual a pseudo-absorption line at 1.785 keV, which is due 
to a calibration problem of the quantum efficiency just 
below the Si absorption edge at 1.839 keV. We suppress 
the corresponding energy channels for the flare phase fit- 
ting, which improves greatly the goodness-of-fit. This cal- 
ibration problem is not visible in spectra with only a few 
thousand counts, and we will neglect it in H5.2.2I The flare 
spectrum is also nearly featureless, with only the He-like 
line of the Fe XXV visible at 6.7 keV. Approximatively 95% 
of the X-ray emission comes from a high energy plasma 
at kT =5.1 keV. We detect a lower energy plasma, which 
is similar to the one found during the quiescent phase, 
with kT - IkeV and EM - 0.2 1 54 c m - 3 . The flare 
abundance is also comparable with the quiescent one; we 
detect no metallicity enhancement as found in the bright 
X-ray flares of active stars (e.g., IGudel et al. 200l|l . We 
obtain during the flare phase an accurate measurement 
of the absorption, A^h = 2.3-3.0 10 21 cm -2 corresponding 
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channel energy (keV) 



Fig. 8. Chandra ACIS-I spectra of LkHa312. The upper 
panel displays separately the ACIS-I data for the quies- 
cent (circle mark) and flare phases. The lines are the cor- 
responding models, a two-temperature plasma combined 
with photoelectric absorption (see Table0]for details). The 
lower panel shows the residual x °f the modeling. 
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Fig. 9. Time-dependent flare spectra. Time interval spec- 
tra are plotted from bottom to top with their one- 
temperature fit (see Table 0] for details) multiplied by an 
arbitrary constant. 



to A v ~1.3-1.7mag UPredehl & Schmitt 1995)l . which is 
consistent with the optical extinction obtain previously 
from SED in 

Thanks to the high count rate of this flare we can make 
a time-dependent spectroscopy of the rise phase to inves- 
tigate the physical properties of this event. 
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Fig. 10. Flare parameters versus time (crosses) compared 
with the quiescent preflare values (grey stripes) from val- 
ues in Table 0] The measured values for the quiescent 
phase are indicated with a dashed-dot line, with a range 
shown as a grey stripe. Time sampling is defined in Fig. EI 
Since the flare spectra can each be fitted with a single 
temperature, we take for comparison with the quiescent 
preflare temperature the average of the two plasma tem- 
peratures weighted over the emission measures, and for 
the quiescent preflare emission measure the sum of the 
two emission measures. 

5.2.2. Time-dependent spectroscopy of the flare phase 

Fig. El shows the flare spectra for the time intervals Ati- 
Atg , which all have a duration of 1 h at the exception of the 
shorter interval Atg having only about 25 min. This bin se- 
lection is a compromise between larger bins, which would 
smooth short term phenomena, and shorter bins, which 
would reduce the number of count. One- hour bins are 
long enough here to collect ^1,500 counts and hence have 
enough statistics to derive accurate physical parameters. 
The spectra are featureless, there is not enough statistics 
to detect the 6.7 keV iron line, which becomes appearent 
only by adding all these time intervals (see Fig. |HJ. We 
find a reliable fit only with a high temperature plasma. 
We tried to add a low temperature plasma having the 
same properties than the one found during the quiescent 
state, but this did not improve significantly the fitting 
result. The X-ray emission of the flare plasma is hence 
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dominated by this high temperature component. However 
as we saw previously in ij5.2.1l bv adding all the flare time 
intervals, it is possible with enough statistics to constrain 
this low temperature emission. We list all the fit parame- 
ters in Table El 

Fig. ^3 plots the flare spectral parameters with their 
errors at the 90% confidence level versus time compared to 
their quiescent values. The absorption is consistent during 
the flare with its quiescent value. The plasma abundance 
is also consistent during the flare with its quiescent value 
with maybe the exception of the interval At3 where the 
metallicity is increased by a factor of 10. This metallicity 
enhancement could be explained by photospheric evapo- 
ration produced by the flare as seen in active stars (e.g., 
Giid el et al. 2 001 ). However assuming that the metallic- 
ity keeps its quiescent value during the flare, this deviant 
bin is consistent with being a statistical fluctuation, as in 
10 independent measurements one discrepant value is ex- 
pected at the 90% confidence level. The emission measure 
increases until At§, an d saturates. 

The emission measure behaviour is hence totally 
different from a rotational modulation of a flare de- 
cay (St elzer et al. 1999jl . where the emission measure de- 
creases exponentially. The plasma temperature is greater 
than the average quiescent temperature during all the 
flare, and reaches its maximum for Ats with kT — 7.6 keV. 
We will consider this maximum as the peak flare in the fol- 
lowing section. The temperature interval Ats-Atg is rem- 
iniscent of the cooling phase of an impulsive flare, which 
displays both an exponential decrease of the temperature 
and the emission measure, however here no decay is visible 
for the emission measure. 

The luminosity in the energy range 0.5-8 keV corrected 
for absorption is Lx.intr ~ 1 X 10 32 ergs -1 throughout 
nearly all the flare. The X-ray luminosity of the flare peak 
reaches ~2% of the stellar bolometric luminosity. The to- 
tal energy released by this flare is Ex — 7.2 1 33 erg in 
8.5 h. 

6. Discussion 

The LkHa312 flare detected here is extraordinary in its 
X-ray luminosity and, due to the serendipitous nature of 
the off-axis location, in the quality of the X-ray data. It 
thus provides a valuable opportunity to examine the ap- 
plicability of established theories for solar-type magnetic 
reconnection flaring in the limit of high power and long 
duration. Note that we consider the absorption, consis- 
tent with a constant value of TVh ~ 1.710 21 cm -2 , to be 
external to the flare environment and of no interest here. 

6.1. Quiescent emission 

The low X-ray emission prior to a major flares, such 
as we see in the first 5 hours of the Chandra obser- 
vation, is usually considered to be a 'quiescent' plasma 
component distributed over large regions as in the solar 
corona. But the quiescent emission here is much hotter 



and more luminous than the solar corona, with the bulk 
of temperatures ranging from 10-35 MK with an aver- 
age (weighted over the emission measure) of 29 MK (see 
H5.2.1I Tabled}. In contrast, the non-flaring solar corona 
has average (weighted over the emission measure) temper- 
ature of ~ 2MK IjPeres et al. 2000|l . The intrinsic quies- 
cent X-ray luminosity in LkHa312 is 6 10 30 erg s _1 in the 
Chandra band (0.5-8 keV) and probably ~ 10 31 ergs -1 in 
a full bolometric X-ray band. In contrast, the bolometric 
X-ray luminosity of the non-flaring solar corona is only 
310 27 ergs -1 during solar minimum to 3 10 28 ergs _1 dur- 
ing solar maximum IjPeres et al. 2 000 ) . Note that the sur- 
face area of LkHa312 is only 6.8 times that of the Sun 
(see H5.2.1(l . so the X-ray luminosity difference of a fac- 
tor of order 1,000 cannot be attributed to a larger surface 
area. 

The higher temperature of the LkHa312 quiescent 
X-ray emission could be attributed to enhanced coronal 
processes. Theories for solar coronal heating based on 
the release of magnetic energy (|Parker 19791 IParker 1 988 
IPriest fc Forbes 2 000 1 have kT ~ fB 2 /(8irn c ), where n c 
is the coronal electron density and / is the fraction of the 
coronal magnetic field which is converted into heat. A 4- 
fold increase on average magnetic field strength, for exam- 
ple, could account for the temperature difference. But it 
seems unlikely that the enormous difference in luminosities 
can be explained in this fashion. Since Lx. quiescent oc n 2 V c 
where V c is the coronal volume, and since n c cannot 
be significantly raised without lowering the temperature, 
a ~ 10 3 increase in coronal volume would have to be 
present to give the 10 31 ergs _1 quiescent luminosity. The 
corona would have to extend >10 stellar radii without sig- 
nificant density decrease, which seems quite unrealistic. 
Alternatively the X-ray volume filling factor in the corona 
of LkHa 312 could be 1,000 times higher than in the solar 
corona. 

The quiescent emission may also arise from the sum 
of many eruptive flares while not excluding a contribution 
from coronal heating processes. This is suggested by the 
growing evidence for microflaring as the origin of most 
of the quiescent X-ray emission from magnetically active 
stars (fOiidel et al. 20031 and references therein). At least 
10 flares each with Lx, flare < 10 30 ergs -1 must continu- 
ously be active over the stellar surface. The lack of vari- 
ation in flux or temperature during the quiescent phase 
suggests a large number of small and moderate-sized flares 
are present, which in turn suggests a dense coverage of 
the stellar surface with high-strength magnetic fields and 
hence a high surface filling factor. This is consistent with 
the fact that the quiescent X-ray luminosity is found to 
reach the saturation level as shown in M5.2.1I 

As the quiescent level appears remarkably stable 
and bright, our understanding of the strong microflar- 
ing of LkHa312 could be advanced by a follow-up X- 
ray observation with high spectral resolution. In par- 
ticular, the density of the flare plasmas can be esti- 
mated using He-like line ratios. NelX at 0.92 keV is sen- 
sitive to densities from a few 10 10 to a few 10 12 cm -3 
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( |Porquet fc Duba"u~ 2000 ) . However, the ionization frac- 
tion of Ne IX may be low due to the high temperatures : 
only 0.4% for plasma at fcT= lkeV and 0.007% at 3keV 
flArnaud fe Rothenflug "1985^ . MgXI at 1.34keV with a 
ionization fraction of 6.8% at fcT=lkeV should then be 
easier to detect, but it probes only higher density from 
a few 10 11 to a few 10 13 cm~ 3 flPorquet fc Dubau 2000| >. 
While the current mission gratings with CCD detectors 
are maybe not sufficiently sensitive for this investigation, 
it will be possible with the forthcoming X-ray spectrome- 
ter aboard ASTRO-E2, a next generation instrument us- 
ing X-ray bolometers. 

6.2. Rise phase 

The morphology of the LkHa 312 flare rise is distinctively 
slow (r ~ 120 min) in contrast to a typical solar flare 
rise (tq ~ 5 min), and smooth with no precursors, shoul- 
ders or secondary events. The smoothness may indicate 
a relatively simple magnetic geometry such as a single 
bipolar loop or bipolar loop arcade without adjacent mul- 
tipolar magnetic areas. The scale length of the loops can 
be tentatively constrained by assuming that the magnetic 
structure confining the X-ray emitting plasma is stable 
during the slow and smooth rise. Loops suffer magne- 
tohydrodynamical instabilities with characteristic growth 
timescale, t, corresponding to the length scale of the pre- 
emption region, I, divided by the ambient Alfven speed 
l|Priest fc Forbes 2000|l . V A oc Bn~ Xj2 , with B the aver- 
age magnetic field strength, and n c the pre-flare coronal 
density. We scale the LkHa 312 values to the solar values : 
; Q = 510 10 cm for an average CME IjGilbert et al. 20011 
and references therein), Bq ~ 100 G (e.g.. IParker 1988|l . 
and n CiQ ~ 10 9 cm -3 (see |[Shibata fc Yokoyama 20021 and 
references therein). We obtain : 

I = l Q (r/r ) (B/B Q ) (n c /n c , e )-^ 2 . (1) 

We point out that this scale length is the scale for the 
volume of erupted field. Generally, it can be larger than 
the volume from which the X-ray emission arises, and in 
the case of the solar flares the scale length of the erupted 
region may be 2 to 3 times greater than the scale of the 
flaring region. Therefore we estimate the following upper 
limit for the height, h, of a semi-circular magnetic loop 
filled by the X-ray emitting plasma : 

h max = 0.5 x 1/3 , (2) 

/i max = 2 10 11 cm x (5/100 G) K/10 9 cm 3 )" 1 / 2 . (3) 

Replacing the loop height by the loop half-length, L = 
irh/2, leads to the following constraint : 

L < 10 8 cm x Bnjf. (4) 

This formula will be used in the next section combined 
with an estimate of L and B, to derive a lower limit esti- 
mate of the pre-flare coronal density. 



6.3. Peak properties 

Both X-ray QFavata fc Micela 2003| and radio 
i|Giidel 2002J1 stellar peak flare properties show that 
magnetic loop scaling relationships can be derived from 
the corresponding solar flare properties. We use here the 
work of Shibata & Yokoyama (1999, 2002) which is based 
on a scaling law relation for the maximum temperature of 
reconnection heated plasma and magnetic field strength 
derived from MHD simulations (Yokoyama & Shibata 
1998, 2001). We note, however, that this work is based 
on MHD simulations where the effect of the radiative 
cooling can be neglected, so scaling laws which include 
the cooling switching from conductive to radiative as 
time proceeds ( |Cargill et al. 1995| ) have not yet been 
numerically tested. 

The theory of Shibata & Yokoyama involves a roughly 
semi-circular magnetic loop with half-length L, magnetic 
field B, and electron density after chromosphcric evapo- 
ration n c . These parameters are derived from the peak 
X-ray temperature and emission measure together with 
prior knowledge of the pre-flare coronal electron density 
n c , using the following scaling formulae where kT is ex- 



pressed in keV : 

L = 2 10 11 cm x EM^l 5 fcT" 8 / 5 n~% 5 , (5) 

B = 32 G x £M 5 - 4 1/5 kT 17 / 10 n 3 J™ , (6) 

n c = 10 10 cm" 3 x EM~ 2/5 kT 12 / 5 n% . (7) 



The loop parameters are hence function of n C) i2, which 
can be constrained. First, for consistency with the chros- 
mospheric evaporation assumption, the flare coronal den- 
sity must be lower than the pre-flare coronal density, which 
leads using @ to the following upper limit estimate for 
n-c,i2 : 

n c ,i2 < 10- 5 EM^-kT 6 . (8) 

Second, from (@J and using _ ® leads to a lower 
limit on n c ,i2 : 

n c ,i2 > 10 9 ^M 5 4 4 fcT- 33/2 . (9) 

As the LkHa 312 peak flare values (EM 54: = 8.0 and kT = 
7.6 keV) are consistent with the empirical correlation (see 
Fig. 1 of |Shibata fc Yokoyama 2002| ), no special scaling 
laws are required for this flare. Applying the previous set 



of formulae, we find : 

0.01 < 710,12 < 0.24, (10) 

L = 2.3 10 10 cm x n~ 2 ^ 2 5 , (11) 

B = 670 Gx n c {™ , (12) 

n c = 6.4 10 11 cm" 3 x n c % . (13) 



Numerically, we find that the loop half-length ranges 
from 4.8-16 10 10 cm (0.3-0.9R*), the loop height ranges 
from 3.1-10 10 10 cm (0.2-0.5R*, i.e., 0.5-1.3 R©), the 
magnetic field from 430-180 G, and the flare loop densities 
from 2.4-0.4xl0 n cm- 3 . 
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6.4. Gradual phase and decay phase 

The gradual phase, which has a slower rise and follows 
the fast rise phase, is really unusual. It may be explained 
by the rapid fall of the stellar magnetic field with height. 
Indeed in this situation when the magnetic loop system 
rises, the reconnection site is continually moving up into 
a weaker field region, which produces the decline of the 
rate of magnetic energy release. 

Much of solar and stellar flare loop modelling is based 
on the decay phase of the flare when the loop plasma is 
cooling (see review bv lReale 20021 e.g.. lReale et al. 1997|) . 
Unfortunately, only the beginning of this phase may have 
been observed in LkHa 312 in the second part of the grad- 
ual phase where we find a decrease of the plasma tempera- 
ture from 88 MK (7.6 keV) to 44 MK (3.8 keV) over 5 hours 
using one-hour averaged spectra (Figure I10|) . However 
during this period, the broad-band luminosity remained 
constant, which is rather unusual for a cooling period, so 
that the emission measure (density and/or volume) in- 
creased somewhat as the plasma cooled, or short heating 
phases may occur which are smoothed by the averaged 
spectra. 



7. Conclusions 

So far, LkHa312 was considered from ~30-yr old Ha 
prism surveys as an unremarkable, weak-activity star at 
the periphery of M 78, until Chandra unveiled its extraor- 
dinary activity in X-rays. Our study based on the available 
optical and IR data shows that it is likely a < 1 Myr-old 
T Tauri star similar to the Sun (M* ~ 0.7-0.75 M , R± ~ 
2.6 R©), about as luminous (Lboi ~ 1-5 L©), and sur- 
rounded by a hollow circumstellar disk. It is therefore 
perhaps a transition object between an active accreting 
'classical' T Tauri star, and a genuine disk-free 'weak-line' 
T Tauri star. 

In contrast with these fairly unoriginal properties of 
LkHa312, our Chandra observations revealed a particu- 
larly long and intense X-ray flare, in addition to a bright 
quiescent emission. Because of the intensity of the flare 
(~ 0.5 Chandra ctss -1 ) and its favorable far off-axis posi- 
tion on the ACIS-I detector which prevented pile-up, the 
quiescent state as well as the rise phase of the flare could 
be studied at an unprecedented level of detail, yielding 
spectroscopy resolved on a time scale of down to 1 h out 
of a 26 h-long observation. The flare itself is quite remark- 
able, with its high X-ray luminosity (comparable to that 
of only a handful of TTS so far) and slow rise, but so is 
the level of the quiescent emission, which is saturated. The 
high intensity of the quiescent level was already percepti- 
ble in the ROSAT/HRI observation of 1997, so it is likely 
not a brief, transient phenomenon. 

Quantitatively, the X-ray properties derived from our 
Chandra observations are as follows : (i) the quiescent 
emission of LkHa 312 is stable over 18 h of observation and 
saturated at a high level of activity (Lx/^bol ~ 10~ 3 , with 
Lx ~ 6 x 10 30 ergs -1 ); (ii) LkHa 312 also displayed an un- 



usually long and intense X-ray flare, showing a long rise 
phase extending over several hours. This flare was bright, 
reaching Lx ~ 10 32 ergs -1 , i.e., ^2% of the bolometric lu- 
minosity. Perhaps more remarkably, the peak temperature 
was quite high (kT ~ 7keV). (iii) In terms of solar-type 
magnetic confinement of the plasma, the application of 
the scaling laws found by Shibata & Yokoyama shows that 
the flare corresponds to a fairly large volume (maximum 
height of about 1 R© ) , and typical solar coronal densities 
(~ 10 n cm- 3 ). 

The problem is to understand the reasons for these 
rather extreme X-ray properties, both in the quiescent 
state and in the flaring state. The first, "standard" way is 
to address the problem in the framework of solar physics. 
This is how we derived the flare properties (size and den- 
sity of the confining magnetic loop, strength of the mag- 
netic field). As discussed in H6.ll one could invoke en- 
hanced coronal processes and/or microflaring, 1,000 times 
more intense than the quiet Sun, or a volume filling fac- 
tor in the LkHa312 corona 1,000 times larger than in 
the solar corona, but the physical reasons for such an en- 
hancement are unclear. This is in fact a general problem 
for pre-main sequence stars : the mechanism for magnetic 
field generation leading to the high level of X-ray activ- 
ity is still not well understood. For instance, the lack of 
an X-ray luminosity/rotation correlation on large sam- 
ples like in Orion is perhaps due to the fact that an or 
(turbulent) dynamo operates at the young TTS stage, 
as opposed to the standard rotation-convection induced 
a— u> dynamo which holds in late- type main-sequence stars 
HFeigelson et al. 20031 1 . 

An alternative way is to consider that the enhanced 
X-ray activity of LkHa312, and perhaps more generally 
of at least some other, very X-ray luminous TTS, is trig- 
gered by an external cause such as a close, faint companion 
like a brown dwarf or even a planet. Indeed, as shown by 
Cuntz et al. (2000)| planets can cause enhanced magnetic 



activity if they are sufficiently massive and close to the 
parent star, like 51 Peg-like planet-bearing stars. The gen- 
eral idea is that tidal effects can enhance the dynamo by 
altering the outer convective zone, and/or, if the planet 
is magnetized, the star and the planet can interact mag- 
netically via shearing and reconnections, for instance as in 
the star-disk magnetic interaction mechanism proposed by 



|Montmerle et al. (2000) to explain the triple X-ray flare of 
the YLW15 protostar. In such a framework, the rotation 
information would be important : although a slow rotation 
would not be conclusive, a fast rotation (periods of days 
or less) could be indicative of a tidally locked, RS CVn- 
like binary system. We note that the circumstellar disk of 
LkHa312 is hollow, with no accretion taking place, which 
is indeed consistent with the possible presence of a planet 
clearing the inner disk fe.g.. INelson 2 003). 

Both ground and space follow-up observations will be 
needed to investigate the coronal properties of LkHa312 
as a genuine TTS and its influence on its circumstellar ma- 
terial : medium resolution optical spectroscopy, to confirm 
its youthfulness with Li I 6707 absorption line measure- 
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ment, to check that the Ha emission has not changed, and 
to estimate its projected equatorial velocity; to determine 
its rotational period from the photometric modulation of 
starspots; high-resolution IR spectroscopy to detect the 
possible presence of a very low mass companion; near-IR 
high resolution imagery to study a possible transition cir- 
cumstellar disk; and of course high resolution X-ray spec- 
troscopy to probe the density of its coronal plasma. 
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